Superconducting thin rings with finite penetration depth 
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Recently Babaei Brojeny and Clem [Phys. Rev. B 68, 174514 (2003)] considered superconducting 
thin-film rings in perpendicular magnetic fields in the ideal Meissner state with negligibly small mag- 
netic penetration depth and presented useful analytical limiting expressions and numerical results 
for the magnetic-field and sheet-current profiles, trapped magnetic fiux, self-inductance, magnetic 
moment, and focusing of magnetic flux into the hole when no net current flows in the ring. The 
present paper generalizes all these results to rings with arbitrary values of the two-dimensional ef- 
fective penetration depth A = A^/d (A is the London depth and d < A/2 the fllm thickness) using a 
straightforward matrix inversion method. We also present results for the energy of a superconduct- 
ing ring as a function of the applied magnetic induction Ba and the quantum number A'' deflning 
the size of the fluxoid N<j)o trapped in the hole. 

PACS numbers: 74.78.-w, 74.25.Ha,74.25.0p 
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I. INTRODUCTION 

In a recent paper Babei Brojeny and Clem- calcu- 
lated the magnetic properties of superconducting thin- 
film rings and disks in the ideally screening Meissner state 
when a perpendicular magnetic field is applied. They 
showed that the Biot-Savart law for the sheet-current 
density J(r) = j'^^^^^'i{r, z)dz (j is the current density) 
in a circular disk of thickness d and radius b with a hole 
of radius a (0 < a < 6) can be solved by a physical 
ansatz containing a factor [(r^ — a^)(6^ — r^)]^^/^ describ- 
ing the divergence of J{r) at the inner and outer radii. 
A plethora of useful numerical results and analytical lim- 
iting expressions is given in Ref. Q for several situations: 
(a) magnetic flux trapped in the hole when no magnetic 
field is applied; (b) zero magnetic flux in the hole when 
the ring is subjected to an applied field; and (c) focusing 
of the magnetic flux into the hole when a magnetic field 
is applied but no net current flows around the ring. 

Throughout the paper^ it was assumed that either the 
London magnetic penetration depth obeys A < d/2 or, 
if A > d/2, the two-dimensional (2D) penetration depth 
A = A^/d (or screening length^ /d) is negligibly small. 
The same assumption was made in previous work on 
disks2i^ and ringsi^ However, while thin superconduct- 
ing rings with A = ideally screen magnetic flux from 
penetrating the hole, a finite A or A = }?/d > A will 
allow magnetic fiux to penetrate into the film as well as 
the hole. This effect is much stronger than would be 
suggested by the exponential factor exp[— (6 — a)/A] that 
applies to long tubes in an axial field. 

The effect of finite A is particularly important for the 
interpretation of experiments that try to confirm the 
cosmological Kibble-Zurek mechanism^'^ of spontaneous 
formation of vortices during rapid cooling of a super- 
fluid; some of these experiments use superconducting 
ringS)^ while others use disksfifi Finite A also modifies 



flux focusing, an important feature of SQUIDs (super- 
conducting quantum interference devices), which usually 
have roughly the shape of a washer with a small central 
hole but with a slit that allows magnetic flux to enter 
the hole and causes zero net circulating current at the 
squid's critical current The concentration of 
magnetic flux into this hole increases the effective area of 
the SQUID. 

This paper is organized as follows. In Sec. II we de- 
scribe our calculation method, which applies to arbitrary 
A. In Sec. Ill we compute the self-inductance of a thin 
flat ring. In Sec. IV we calculate the response of a ring 
in an applied magnetic held. In Sec. V we study the 
flux-focusing problem and calculate the effective area. In 
Sec. VI we calculate the energy of a ring as a function 
of the applied field and the quantum number N describ- 
ing the size of the fiuxoid N(j)o trapped in the ring. In 
Sec. VII we give some analytical results for the limit of 
large A 6, which applies to mesoscopic rings. We 
present a brief summary in Sec. VIII. 



II. CALCULATION METHOD 

We assume for simplicity that current in a circular coil 
far from the ring produces a vector potential AQ(r) — 
Aa(v)ip, which describes the magnetic induction Bq = 
V X Aa', If is the azimuthal unit vector. Near the film 
in the plane z = 0, we assume Aa{r) — rBa/2, such that 
the magnetic induction applied to the thin-film ring is 
Bq = BaZ. In response to either the applied field or 
a fluxoid trapped in the hole, currents are induced in 
the ring. The net magnetic induction is B(r) = BQ(r) + 
Bj(r), where Bj(r) = VxAj(r) and its vector potential 
Aj(r) are generated by the currents in the film. Because 
of the circular symmetry, the sheet current in the film 
has only a (^-component, 3{x,y) = J{r)ip. Similarly, the 
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vector potential A{x,y,z), defining the total magnetic 
induction B = V x A, has only a yj-component. In the 
film plane z = 0, we have A(a::, y, 0) — A(r)ip and A{r) = 
^air) + Aj{r). The current density j is related to A via 
the London equation, j — — A^Z/igA^, where A^ is the 
superfluid velocity expressed in units of vector potential; 
here A^ = Asip with 



As{r) = A{r) - <I>//27rr. 



(1) 



The second term on the right-hand side is due to the 
gradient of the phase of the complex superconducting 
order parameter. For a ring with a slit we treat $y as a 
free parameter to be determined by boundary conditions, 
but for an unslitted ring corresponds to the London 
fluxoidf^ which is quantized; i.e., $/ = NcfiQ, where N 
is an integer and 00 = h/2e is the superconducting flux 
quantum. 

It is useful to divide the fields into two contributions, 
B = Bi +B2, ^ = + A2, j = ji +j2, J = Ji + J2, 
etc., where the subscript n = 1 indicates that it is 
driven by the fluxoid [driving term Di{r) = —^fj^nr 
and A\{r) = Aji(r)], and the subscript n = 2 indi- 
cates that it is driven by the applied field [driving term 
D2{r) = Aa{r) and A2{r) ^ Aa{r) + Aj2{r)]. The Lon- 
don equation for contribution n is 



J„(r) = -[Dr,{r) + A,j,,{r)]/ iioK. 



(2) 



In this paper we will calculate the two contributions (n 
= 1 and 2) to the magnetic flux $„(a) through the hole, 
total flux through the ring, magnetic moment m^, 

and total current /„ around the ring using the definitions 

$„(r) = 2^ / dr' r'Bnir') = 2^rA„(r) , (3) 
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m„ = TT / dr r J„(r), /„=/drJ„(r). (4) 

J a J a 

The relation between the sheet current J„(r) and the 
vector potential Aj„(r) it generates is obtained as fol- 
lows. From the Maxwell equation /xoj = V x B = 
V X V X A = — V^A (since here V • A = 0) we obtain the 
3D Biot-Savart law for the current-generated part A,7(r) 
of A: 



A,/(r) = ^0 



d^r 



j(r') 



Anlr — r' 



(5) 



with r = {x,y,z). Integrating this over z' and ip, noting 
that j flows only inside the film —d/2 < z < d/2, one 
obtains in the plane z — for each component: 



^Jn(r) -/io / dr'J„(r')Q(r,r') 



with the integral kernel 
r'^dp 



Q{r,r') 



r cos ip 



27r (r-^ + — 2rr' cos (p)!/^ 



(6) 



(7) 



This kernel may be written in terms of elliptic inte- 
grals, but for transparency we prefer a fast direct nu- 
merical integration of (7); see also Refs. 2 3 4 d High 
accuracy is achieved by substituting in the integral p = 
p{u) — TTM — sin(7ru) and integrating over < u < 1 
using an equidistant grid for u, Ui = {i — l/2)/A^^, 
i = 1,2, . . . , N^p, w 30 — 60, with weights Wi — 
ip'{u)/N^ = [1 - cos(7rui)]7r/iVy: 



Q(r,r')= / dpf{p>)^ / duf[p{u)]p'{u) 



Ui Wi 



(8) 



i=l 



Writing J„(r) ~ J dr' Jn{r')S{r — r') and inserting Eq. 
(2) into Eq. (6), we obtain 

A.(r) = -/io / dr' J4r')[Qir,r')+ASir-r')]. (9) 



(For introduction of finite A into other geometries see 
Ref. Formally, the integral equation (9) may be 

solved for the sheet current J{r) by writing 

nb 

Jn{r) = -^lo^ / dr'D„{r')K{r,r') , (10) 



where K{r, r') is the inverse of the kernel Q{r, r')+AS{r- 
r'), defined by 



dr'K{r, r')[Q{r', r") + K5{r' -r")\ = 5{r-r") . (11) 



The inverse kernel K{r, r') is easily calculated numeri- 
cally by introducing an appropriate grid r.^ with weights 
Wi such that the integral is approximated by a sum. 



.6 N,. 
/ dr/(r)«^/(r,) 



Wi . 



(12) 



High accuracy is achieved in the present case, where the 
integrated function /(r) may have infinities at r = a 
and r = 6, by a grid that is very dense near r = a and 
r = b. A good such grid is found by the substitution 
r = r{u) = a+{b-a)il0u^-15u^ + 6u^), r'{u) =30(6- 
a){u — u^Y , yielding = ?'(ui), Wi = r'{ui)/Nr, with 
M, 1/2) /Nr, i^l,2,...,Nr, Nr^SO- 100. This 

grid defines the vectors Jni = Jn(fi), Dni = Dn{ri), and 
the matrix Qij — Q{ri,rj). Equation (9) then becomes 
a sum (or matrix multiplication): 



Nr 



Dm = - ^{wiQij + ASij) ^j,oJ„ 



This is inverted by 



Nr 



3=1 



(13) 



(14) 
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where Kij is an inverse matrix: 

K,j = {w^Q.j + AS,j)-^ (15) 

(no summation over i; dij = I if i — j; otherwise Sij = 0). 
The matrix equation from which the total sheet current 
is determined for given Ba and $/ thus reads expHcitly: 

,oJ(..) = -E^.(^^-2^)- (16) 
j=i J 

Finally, one hint is required without which this method 
may not work or is inaccurate. The matrix Qij has infi- 
nite diagonal terms Qa, since Q(r,r'), Eq. (7), diverges 
logarithmically when r approaches r'. Namely, one has 

1 ri 

Qi,i+i ^ Q{ri,ri +Wi) ^ —In — . (17) 

This problem was dealt with in detail in Refs. [3 and 
U where the optimum choice of the Qu was found first 
numericallj*^ and then analyticallj*^, e.g., from the con- 
dition that an infinite disk ideally screens two coils sepa- 
rated by the disk. One thus has the complete definition 
of the matrix Qij : 

Qu = —(in ^-2 . 18 

In Eq. (15) the diagonal term is WiQa + A. Thus, when A 
is larger than the maximum value of the weight Wi (or the 
spacing between grid points), which is of order (b—a) /Nr, 
then the choice of the Qu is not critical, and for the 
computation of J(r) one may even put Qu = 0. For 
small A, however, the correct choice of Qu is important. 

III. SELF-INDUCTANCE OF FLAT RINGS 

According to Eq. (16) the current in a ring originates 
from either an applied field Ba, trapped flux related to 
the parameter <&/ , or both. In this section we put Ba = 
and compute the sheet current Ji(r) for finite trapped 
fluxoid The total current Ii is then given by Eq. (4). 
The current generates a flux <f>i(r) = 2TTrAi{r), Eq. (3), 
with Ai{r) — Aji{r) from Eq. (6). Figure 1 shows the 
profiles of the sheet current Ji (r) , perpendicular induc- 
tion Bi{r) = {l/r){rAiy , and vector potential Ai{r) for 
a ring with a/b = 0.4 for three A/6 values (0, 0.03, 0.1) 
for the case of flux trapping with Ba = and > 0. 
One can see that finite A removes the infinities of Ji(r) 
at r = a and r = b [Ji oz (r — a)~^/^ and Ji cx (6 — r)~^/^; 
see Ref. which lead to a similar infinity of Bi{r). For 
A > 0, Ji(a) and Ji{b) are finite, and Bi{r) penetrates 
the superconductor but still exhibits a logarithmic infin- 
ity at r = a and r = b, which is caused by the abrupt 
jump of Ji{r) to zero. Note also that with increasing A 
the total current Ii, the magnetic flux ^i{a) in the hole, 



3 



2 



< 

CM 

-1 
CD 







-1 





^1 ' 

/ ] 


1 \ ' 

\ Vi 






" A/b = . 
0.03 
0.1 ^ 


/ 1 

\/ ' 
/\ 1 \ 

X ' 


1 \\ 






















3 0.2 


0.4 0.6 


0.8 


1.2 



r/b 



FIG. 1: Profiles of sheet current Ji(r) in the ring, perpen- 
dicular magnetic induction Bi{r) = (l/r)(r^i)', and current- 
generated vector potential Ai{r) = Aji{r) in the plane of the 
ring for the case of trapped flux and zero applied field {Ba = 0, 
>!>/ > 0, Figs. 1-6). The hole radius is a = OAb and the 2D 
penetration depth is A = (solid curves), A — O.OSfe (dot- 
dashed), and A — 0.1b (dashed). Plotted are the dimension- 
less quantities Ji/($///io&^), -Bi/(<E>//6^), and 2Ai/{^f/b). 
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FIG. 2: Computed self-inductance L = <l"///i of a supercon- 
ducting thin flat ring with hole radius a and outer radius b, 
plotted in units of noa for 2D penetration depths A — 0.001, 
0.01, 0.03, 0.1, and 0.3 in units of 6. Note that even very 
small A <^ b considerably enhances the inductance of rings 
with a small central hole {a/b ^ 1) as compared with the 
case A = 0. The circles depict the empirical expression (20) 
valid for A = 0; they perfectly coincide with the dots that 
mark the curve computed for A = 0. The dashed curve shows 
the limiting expression Li, Eq. (19). See also Fig. 6 below, 
showing the inverse self-inductance fioa/L. 



FIG. 3: Magnetic flux "l>i (a) in the hole of a superconducting 
thin flat ring with hole radius a and outer radius b for the 
trapped-flux case as in Figs. 1-6 (applied field Ba = and 
fiuxoid <!>/ > 0) for 2D penetration depths A/b = 0, 0.0001, 
0.001, 0.01, 0.03, 0.1, 0.3, 1, 3, and 10. The dots on some of 
the curves mark the a/b grid used here in all such figures. 



FIG. 4: Minimum of the magnetic field -Bi(O) occurring in 
the center of the ring for the trapped-flux case {Ba = and 
$/ > 0) plotted as 27ra^Bi(0)/$/ for A/& = 0, 0.0001, 0.0003, 
0.001, 0.003, 0.01, 0.03, 0.1, 0.3, 1, 3, and 10. The dashed 
curves show the large-A approximation, Eq. (41), for A/b — 1, 
3, 10. 



and the field Bi in the hole decrease^ while the London 
fliDcoid 2Trr[Ai{r) + /ioAJi(r)] = <&/ remains constantiii 
In the case of ideal screening, i.e. for A = 0, no mag- 
netic field penetrates the ring material, and thus the flux 
$i(r) is the same for any r between r — a and r — b. 
The self-inductance of the ring may then be defined as 
L = <i>i//i with $1 = $i(a) = $1(6). For small reduced 
inner radius a = a/b 1 the flux trapped in the hole was 
found by Clemiiii to be <i>i(a) — thus the induc- 

tance L = ^i/Ii approaches Lq = 2/ioa- In the opposite 
limit, for narrow rings with a ^ 1, L approachesii^ 

Li = ^oi?[ln(8i?/u;) -2 + ln4] 

= ^o&(l + a)(tanh"^a- l + ln4) (19) 

[R = {a + b)/2, w = b — a]. For arbitrary a, but still 
A = 0, L was computed in Ref. 1, where a useful empirical 
formula was presented, 

L2 = fiob[a -0.197a'^-0.03la^ + {l + a)t'Anh^^a] . (20) 

This fit is confirmed by our method, as shown in Fig. 2; 
its relative deviation from the exact L ranges from 
-0.005% to -HO.06%. 

When the effective penetration depth A is finite, the 
magnetic field penetrates into the ring material, and the 
magnetic flux $i(r) is no longer constant when r changes 
from r = a to r — b. The definition of L = ^i/Ii 
via magnetic flux is thus not unique; in particular, 
$i(a)//i (f>i(&)//i when A > 0. In this general case 
one must use the definition of L via the electromagnetic 
energy of the ring, E — {l/2)Llf. This energy is com- 
posed of the magnetic energy Em and the kinetic energy 
Ek of the currents, as is evident from the energy integral 




FIG. 5: Magnetic moment mi of a superconducting thin flat 
ring with inner and outer radii a and 6 for the trapped-flux 
case {Ba = and $/ > 0) plotted as the dimensionless ratio 
= nomi/b^f for A/b = 0, 0.001, 0.01, 0.03, 0.1, 0.3, 1, 3, 
and 10. The dashed curves show the large-A approximation, 
Eq. (42), for A/6 = 1, 3, 10. See also Fig. 12 below. 

of London theory: 

E = Em + Ek = ^ j(fr{nl+XY,) (21) 

with Hi = ytio ^Bi and ji = V x Hi. The integral (21) 
over all space can be transformed into an integral over 
the superconductor by introducing the vector potential. 
For a flat ring this yields 

E^Em+Ek^TT f drr[Ji(r)Ai(r) +//oAJi(r)2].(22) 

J a 



FIG. 6: Total current Ii in the flat ring with trapped flux 
{Ba = and $/ > 0) plotted as the dimensionless ratio aj = 
Hoah/^f = noa/L for A/6 = 0, 0.0001, 0.001, 0.01, 0.03, 0.1, 
0.3, 1, 3, and 10. Same data as in Fig. 2 but inverted. The 
dashed curves show the large-A approximation, Eq. (40), for 
A/6= 1, 3, 10. 

From the two energy terms when Ba = one may de- 
fine the geometric inductance L,„ = 2Em/Ii and the 
kinetic inductance Lk = 2Ek/ll, yielding the total self- 
inductance L = Lm + Lk = 2E/If. The energy E, 
Eq. (22), may also be written as 

E = E^ + Ek = il/2)<Pfh, (23) 

where $/ = 2Trr{Ai + /zqAJi) — const (for a < r < b) 
is the London fluxoidii^iii We may derive Eq. (23) from 
(22) by noting that the constant combination 27Tr{Ai + 
fj,oAJi) = can be factored out of the integrand. 

From Eq. (23) and the definition E = {l/2)Llf it im- 
mediately follows that 

L = (24) 

This general result, valid for arbitrary penetration depth 
A, differs from the definition used previously^ for A = 
0, in that the magnetic flux <i>i(a) through the hole is 
replaced by the fluxoid $ / , which coincides with the flux 
<i>i(a) in the hole only in the special case A = 0. 

Equation (23) also may be derived by considering how 
the fluxoid in the ring may be increased from zero to $/ 
by moving vortices (Pearl vorticesi2ii242& of short length 
d) from the outer radius through the ring into the hole. 
Each vortex has to cross the current-carrying ring, where 
a Lorentz force i^q Ji(r) acts on it. Integrating this force 
from r = b to r = a, one obtains the energy ^q/i. Each 
crossing vortex increases the phase change of the super- 
conductor order parameter around a circle in the ring by 
27r and thus increases the fluxoid by (pa . Noting that the 
total current Ii is proportional to the number of vortices 
that already are in the hole, one obtains E = (l/2)$y/i, 
Eq. (23). 



FIG. 7: Profiles of sheet current J2{r) in the ring, perpen- 
dicular magnetic induction B2{r) — (l/r)(rA2)', and vector 
potential A2{r) in the plane of the ring for the zero-fluxoid 
state in a finite applied field {Ba > 0, $/ = 0, Figs. 7-12). 
The hole radius is a = 0.46 and the 2D penetration depth is 
A = (solid curves), A = 0.036 (dot-dashed), or A = 0.16 
(dashed), as in Fig. 1. Plotted are the dimensionless quanti- 
ties J2/Ha, B2/Ba, and ~8A2/bBa. 

Figure 2 shows the inductance L — ^f/Ii for various 
ratios A/6 = 0, 0.001, 0.01, 0.03, 0.1, and 0.3. The circles 
show the fit (20), which is an excellent approximation 
for A = and all hole radii a. Note that L increases 
with increasing penetration depth A. Even small A/b — 
0.001 noticeably enhances L of rings with a small hole. 
Figure 3 shows the magnetic flux in the hole, <I>i(a) = 
2iTaAi(a) , referred to the fluxoid <&/ trapped in the ring. 
Note that $i(a) may be much smaller than <i>/ even for 
small A/b. A similar plot. Fig. 4, shows the minimum 
field in the hole, Bi{0), occurring at the center r — 
(cf. Fig. 1) plotted as 27ra^Bi(0)/$/. These curves look 
qualitatively similar to <i>i(a)/$/ in Fig. 3. For A = 
and a ^ & one has 27ra^i?i(0)/<i>/ — > 1. 

The magnetic moment mi of this ring (still for Ba = 
and $/ > 0) is depicted in Fig. 5 as the dimensionless 
ratio am — fJ-o'nT'i/b^ for further approximations to 
these curves see the dotted and dot-dashed curves in the 
similar Fig. 12 below. The total current Ii times the 
inner radius a is depicted in Fig. 6 as the dimensionless 
ratio a J = fiQali/^f. Actually Fig. 6 shows the same 
data as Fig. 2, but inverted, since the plotted quantity 
liQali/^f equals fiaa/L; however. Fig. 6 shows the entire 
range < A < oo. 



IV. ZERO-FLUXOID STATE 

This section considers the zcro-fluxoid state reached 
when a thin superconducting ring is slowly cooled in zero 
field. A magnetic field Ba = fioHa is then applied, but 
the fiuxoid ^ f in the ring then remains zero. When A = 
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FIG. 8: Magnetic flux $2(0) in the hole of a flat ring when 
Ba > and = in units of -Ka^Ba for A/6 = 0, 0.0001, 
0.0003, 0.001, 0.003, 0.01, 0.03, 0.1, 0.3, 1, 3, and 10. The 
bottom plot shows the same data ten times enlarged. 



0, one has ideal screening, such that the magnetic flux 
$(a) in the hole is zero. For flnite A > the screening 
is incomplete, some flux leaks through the ring, and the 
flux $(a) in the hole is no longer zero. In this section 
we therefore put Ba > and <&/ = in our numerics. 
Figure 7 shows for this case the proflles of J2{r), B^ir), 
and A2{r) for a ring with a/b = 0.4 for K/h = 0, 0.03, and 
0.1. Note that the induction B2{r) changes sign inside 
the hole and has a negative infinity at the inner edge 
r = a and a positive infinity at the outer edge r = b. 
For A = the integral of B2{r) over the hole area is 
$2(0) = 0, but for A > the flux $2(0) > 0. 

Figure 8 shows the penetrated flux $2(0) in units of 
its maximum value ira^Ba reached in the limit of A ^ 6, 
and Fig. 9 shows the fleld maximum i?2(0) occurring at 
the center r = of the ring. Both $2(0) and i?2(0) in- 
crease monotonically with both a and A, but while ^2(0) 
at small a increases linearly with the radius a, the pen- 
etrated flux $2(0) at small a and small A has negative 
curvature; see bottom of Fig. 8. The penetrated flux 
$2(a)/7ra^Ba ~ B2{Q)/Ba in the limit of small hole ra- 



FIG. 9: Magnetic fleld £2(0) at the center r = of a thin 
flat ring when Ba > and $/ = in units Ba for A/b = 0, 
0.0003, 0.001, 0.003, 0.01, 0.03, 0.1, 0.3, 1, 3, and 10. The 
dashed curves show the large- A approximation, Eq. (43), for 
A/b = 1, 3, 10. 
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FIG. 10: Penetrated flux in thin flat rings with small hole 
radius a <^ b for Ba > 0, $/ =0, plotted versus A/b as 
^2(0) /na^Ba (dots), which in this limit is nearly equal to 
B2{0)/Ba (solid curve) and is well approximated by Eq. (25), 
shown as the dashed curve. 



dius a/b <C 1 is depicted in Fig. 10 as a function of A/b. 
This curve is well fitted by 



$2 (a) B2{0) 



i-hitanh(^2.881n^ -0.675) . (25) 



ira^Ba Ba 2 2 



Figure 11 shows the magnetic moment m2 of the ring 
in units of mg — —{8/3)b'^Ha, which is reached for ideally 
screening disks (a = A = 0);^'^^ the dimensionless ratio 
plotted is (3rn = m/niQ. For A = and a ^ b one has the 
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FIG. 11: Magnetic moment 7712 of a thin fiat ring when 
Ba > and $/ — 0, plotted as the dimensionless ratio 
f3m = m2/mo, where mo = -{8/3)Hab^ for A/6 = 0, 0.001, 
0.01, 0.03, 0.1, 0.3, 1, 3, and 10. The dot-dashed curve shows 
the limit 6 — > o, Eq. (26). The dashed curves show the large-A 
approximation, Eq. (44), for A/6 = 1, 3, 10. 




a/b 



FIG. 12: Total current I2 in a thin flat ring when Ba > 
and $/ = for A/6 = 0, 0.001, 0.003, 0.01, 0.03, 0.1, 0.3, 1, 
3, and 10, plotted as the dimensionless ratio fii = — /2/6-ffa. 
The dotted and dot-dashed curves show the limits, Eq. (27). 
Note that the curves of /3i — —li/bHa agree exactly with 
those in Fig. 5, am = fiomi/b^f , although they describe 
different physical quantities. The dashed curves show the 
large-A approximation, Eq. (45), for A/6 = 1, 3, 10. 
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{l + a/bf 



128tanh-^(a/6) - 1 -I- In 4 



(26) 



Figure 12 shows the total current I2 in the ring induced 
by the applied field, expressed in terms of the dimension- 
less ratio Pj = —l2/bHa. The limits for A = are^ 
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FIG. 13: Profiles of sheet current J(r) in the ring, perpen- 
dicular magnetic induction B{r) = (l/r)(ryl)'m and vector 
potential A{r) in the plane of the ring for the case of zero to- 
tal current / = (flux focusing. Figs. 13-17). The hole radius 
is a = 0.46 and the 2D penetration depths are A = (solid 
curves), A = 0.036 (dot-dashed), or A = 0.16 (dashed), as in 
Figs. 1 and 7. Plotted are the dimensionless quantities J/ Ha, 
B/Ba, and ~4A/bBa. 



(tt/A) (a + b) Hg 
tanh"\a/5) - 1 -fin 4 



for 



(27) 



-{4/7T)bHa 



for a ^ , 



Note that the curves in Fig. 12 exactly coincide with the 
curves in Fig. 5, even though they depict different physi- 
cal quantities for different cases. This identity, = /3/, 
can be proved by evaluating the sum of energy integrals 
P12 = /(i'^rBi-B2//io+/c''^f/^oA^ji-j2 , where the first in- 
tegral extends over all space, including the coil producing 
the applied field Ba, and the second integral extends only 
over the volume of the ring. With the help of the vector 
identity V-(AxB) = B-Vx A-A-(VxB) with A = A2 
and B = Bi, the divergence theorem, and Eq. (2) with 
n = 2, we can show that i^i2 ~ 0. Then, using the same 
vector identity but with A = Ai and B = B2, the diver- 
gence theorem, Ampere's law, and Eqs. (2), (4), (6), and 
(7) with n = 1, we obtain F12 — miBa + ^ fl2 = 0, which 
yields a™ = Mow 1/6$/ — —fiQl2/bBa — Pi- This also 
can be proved by inserting Ji and J2 of Eq. (10) into the 
definitions (4) of toi and I2, renaming the variables r <-> 
r' , and noting the symmetry of K(r,r') — K{r',r)(r' /r) 
defined by Eq. (11). This symmetry follows from the 
symmetry of Q{r,r') = Q{r',r){r'/r) defined by Eq. (7), 
from Eq. (11) noting that [r"lr)5{r - r") 5{r - r"), 
and from the additional property of the inverse kernel 
that Sldr'[Q{r",r') + K5{r" - r')]K{r',r) = 5{r" - r). 



V. FLUX FOCUSING 

This section considers a ring that is in a perpendicular 
magnetic field Ba and contains a fiuxoid <i> f chosen such 
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FIG. 14: Reciprocal of the effective area AcS = ^f/Ba of a 
thin fiat ring with / = (flux focusing) plotted as (nab) /Aea 
versus the radius ratio a/b for penetration depths A/6 = 0, 
0.0001, 0.0003, 0.001, 0.003, 0.01, 0.03, 0.1, 0.3, 1, 3, and 
10. The circles mark the parabolic fit, Eq. (29), and the dots 
nearly coinciding with the curves A/6 = 3 and 10 show the 
limit A ^ oo, Eqs. (30), (47). 




FIG. 15: Magnetic flux ^(a) in the hole of a thin flat ring 
with / = (flux focusing) plotted as ^{a)/iTabBa versus a/b 
for A/6 = 0, 0.0001, 0.0003, 0.001, 0.003, 0.01, 0.03, 0.1, 0.3, 
1, 3, and 10. 



0.2 0.4 , , 0.6 0.8 1 

a/b 

FIG. 16: Magnetic induction B{0) in the center of a thin 
flat ring with / = (flux focusing) plotted as [B{0)/ Ba]{a/b) 
versus a/b for A/6 = 0, 0.0003, 0.001, 0.003, 0.01, 0.03, 0.1, 
0.3, 1, 3, and 10. The dashed curves show the large-A ap- 
proximation, Eq. (48), for A/6 = 0.3, 1, 3, 10. 




FIG. 17: Magnetic moment m of a thin flat ring when I — 
(flux focusing) in units of mo = — (8/3)ffa6"^ for A/6 = 0, 
0.0003, 0.001, 0.003, 0.01, 0.03, 0.1, 0.3, 1, 3, and 10. The 
dashed curves show the large-A approximation, Eq. (49), for 
A/6 = 0.3, 1, 3, 10. 





that no net current circulates around the ring. This cir- 
cularly symmetric situation approximates a ring with a 
narrow slit along one of its radii as proposed by Clemii*ii 
It neglects both the magnetic flux in the slit and the ef- 
fects of the radial currents that flow in opposite direc- 
tions alongside the slit. The slit interrupts the current I 
and allows magnetic flux to penetrate into the hole such 
that the condition / = holds. When Ba is increased 
with ramp rate dBa/dt, a voltage U appears across the 
slit at radius i, U — —d^j/dt — —{dBa/dt)Aeg, where 
^cff — /Ba is the effective area of the slitted ring. In 



our approximation of a circularly symmetric ring the flux- 
oid <&/ is given by = 2Trr[{r/2)Ba + Aj{r) + fj,oAJ{r) ] 
for any r in the superconductor, a < r < b, and it equals 
the parameter <I>^ entering Eq. (2). In a superconductor 
without a slit, <&/ = Nc/iQ is quantized, with integer N. 
When the ring has a radial slit, it may be used to con- 
struct a dc SQUID by connecting the two banks of the 
slit to a superconducting current lead via two identical 
small Josephson junctions, each with maximum super- 
current IgiiSii^ The critical current of the resulting dc 
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SQUID is 

h - 2/0 1 cos(7rB,^ff/0o)| , (28) 

where the effective area is ^cff — ^f/Ba- Note that the 
effective area can be calculated as Acs = ^{a)/Ba (as in 
Ref. Q) only in the limit A = 0, when the magnetic flux 
$(a) in the hole of the ring is exactly equal to the fluxoid 
$ / . The term "flux focusing" is appropriate because Agg 
is always larger than the actual area of the hole (7ra^), 
regardless of the value of A; moreover, Aes is always in 
the range tto^ < AcS < irb^. 

When A = 0, one has the limits^ Aoff/(7ra^) = 
(8/7r2)(6/a) for a < & and Aes/iTra^) = l + {l-a/b) for 
a ^ b. A good flt valid for all a in the range < a < 6 
is (see Fig. 14) 

Aeff/(™2) = l/M 1 + (vrVS - 1)(1 - a)']} (29) 

with d = a/b. Values of Agg obtained from this ex- 
pression deviate bv less than 0.5% from the A = cal- 
culations of Ref. lili^ In the limit of large A ^ b (which 
may be applicable to mesoscopic ringa^iSi^) one has (see 
Eq. (46) below) 

Aeff/(7ra2) = {b^ ^a^)/[2a^ln{b/a)], (30) 

which is already closely approached for A/6 > 1 (see 
Fig. 14). 

The condition that I = h + I2 =0 (or /2 = -/i) 
yields Acs = abj3i/ai, where /?/ = —fiol2/bBa, shown in 
Fig. 12, was computed in Sec. IV and aj = fioali/^f, 
shown in Fig. 6, was computed in Sec. III. Figure 13 
shows the profllcs of J(r), B{r), and A{r) for the case of 
flux focusing in a ring with a/b = 0.4 for A/b — 0, 0.03, 
and 0.1. Note that the sheet current J(r) changes sign 
inside the superconductor and has zero integral, / = 0. 
As in Figs. 1 and 7, when A = 0, J(r) has inverse square- 
root infinities at r = a and r = 6, which were treated in 
Ref. but for A > 0, J(a) and J{b) have finite values. 

In Fig. 14 the reciprocal of the effective area Acs is 
plotted in the form nab/Acs, such that the data for all 
values of the penetration depth < A/6 < cxo can be 
presented in one plot. The curve for A = is well fitted 
by a parabola, Eq. (29), ranging from tt^/S at a/6 ^ to 
1 at a/b —^ 1. The limiting curve for A/6 > 1, Eq. (30), 
is practically reached already when A/b exceeds unity. 
The main message from this plot is that A^e increases 
with increasing A for any given a/6, and that for A 7^ 0, 
Acs/nab diverges when a/6 0, while for A = it tends 
to the finite value S/tt^. For A/6 < 0.03, Af.^/irab has a 
minimum (since 'Kab/Acs has a maximum) as a function 
of a/6. 

Figures 15, 16, and 17 show the magnetic flux ^(a) in 
the hole of the ring plotted as ^{a) / nabBa, the mag- 
netic field B(0) in the center of the ring plotted as 
[B{Q) / Ba\{a/b), and the magnetic moment m/niQ of the 
ring for the fiux-focusing case / = 0. All three quantities 
decrease with increasing A/6. 



VI. RING ENERGIES 

The energy of a superconducting ring with A ^ 6 was 
calculated in Ref. 24 as a function of the fluxoid number 
N and the applied magnetic induction Ba- Using the 
results of the previous sections, we are able to calculate 
the energy of a thin ring {d < A) with inner and outer 
radii a and 6 for any value of A. The authors of Ref. 
also calculated, assuming A 6, the energy barriers 
between states N and iV ± 1 associated with the energy 
cost of moving a vortex or antivortex between the inner 
and outer radii. However, to extend such calculations to 
the case of arbitrary A is beyond the scope of our paper. 

We begin by calculating the total electromagnetic en- 
ergy of the ring-coil system, where the coil produces a 
perpendicular magnetic induction Ba at the ring, 

E = J d^r B^/2nn + y"dV^oA^jV2. (31) 

The hole contains the fluxoid = N(j)o, and B = 
Bq-I-Bj, as discussed in Sec. II. The first integral, which 
extends over all space, is the total magnetic-field energy, 
and the second integral, which extends only over the vol- 
ume of the ring, is the total kinetic energy of the su- 
percurrents. Using the vector potential, the divergence 
theorem, and Eqs. (2) and (4), we obtain 

E^Ea+ mBj2 + N(j,aI/2, (32) 

where Ea = J d^r B^/2/io is the applied field's magnetic 
energy in the absence of the ring, m = mi + ni2 is the 
total magnetic moment of the ring, and I = I1+ 12 is the 
total current around the ring. When the coil currents are 
controlled so as to maintain the magnetic induction Ba 
applied to the ring, the relevant energy is the Gibbs free 
energy, which we define as GN{Ba) = Ea + gN(,Ba) = 
E — W, where W = rriBa is the work done by the power 
supply in bringing up the magnetic moment to its final 
value m. Thus 

9N{Ba) = -mBa/2 + N(f>nI/2. (33) 

Eliminating m and / in favor of the quantities we have 
defined and calculated in Sees. II-V (mi = 6iV(/)oam//io, 
m2 = -{8/3)b^HaPm, h = N(t>oai/^ioa, h = -&i?aA, 
L = ^Qa/aj, and /?/ = am), we may express gN^Ba) in 
a form equivalent to that given in Ref. M for the case m 
which there are no vortices in the annular region a < r < 
6: 

gNiBa)^eo[ih~N)^+jh'']. (34) 

The characteristic energy is eq = 0o/2-^ = 0Qa//2/zoa. 
As can be seen most clearly from Fig. 2, for given val- 
ues of a and 6, eq decreases monotonically as A increases 
from zero to 00 (see Sec. VII); accordingly, eq is a mono- 
tonically decreasing function of temperature. 

The reduced field in Eq. (34) is h = Ba/Bo, where the 
scaling field is Bq = ipo/Acs, with A^g = ab/Sj/aj. Since 
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FIG. 18: Coefficient 7 of ft^ = {Ba/Bof in the A^- 
independent part of the energy of a superconducting ring in 
an applied field [see Eq. (34)], plotted as 1/(7 + 1) = 1/x to 
show the entire range of < a/b < 1. The dots mark the 
limit A ^ 00, Eq. (51). 
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FIG. 19: A''-dependent part of the energy of a superconduct- 
ing ring in an applied field [see Eq. (34)], plotted as (h — N)"^ 
vs the reduced field h = Ba/Bo. 



Bq is proportional to the reciprocal of Aes, shown in Fig. 
14, we see that Bq decreases from its largest value when 
A = to its smallest value when A = cxd. While Bq is 
a monotonically decreasing function of temperature, the 
range of values spanned by Bq is large only for small a/b; 
in a narrow ring for which a b, we find Bq (pQ/nb'^, 
independent of A. 

The constant 7 in the second term of Eq. (34), ob- 
tained as 7 = X ~ 1 = 86/3ma//3aa^ — 1, is shown 
in Fig. 18. This A^-independent quadratic term is rel- 
atively unimportant, however, because it is only the first 
term on the right-hand side of Eq. (34) that determines 
which quantum state N has the lowest energy. Figure 
19 exhibits a plot of {h — iV)^ vs h = Ba/Bo for sev- 
eral values of N. From this plot we can see clearly that 
the state N has the lowest energy for values of h in the 
range N - 1/2 < h < N + 1/2. When h ^ N - 1/2, 
the states N and — 1 have the same energy; similarly, 
when h = A^-|- 1/2, the states TV and iV-|- 1 have the same 
energy. 

The temperature dependence of Bq = <^o/^cff leads 
to the possibility that the energetically favored value of 
N may be higher close to Tc, where A diverges, than 
at lower temperatures, where A is much smaller. This 
effect, which may have important consequences for ex- 
periments using small superconducting rings, is greatest 
when a/b is small. For example, suppose that a — 0.1b. 
At high temperatures for which A/b > 1, we see from 
Fig. 14 that A^g « 0.2nb'^, such that Bq « S^o/tt^^. 
In an applied field of Ba ~ 30o/7r6^, h « 0.6, and the 
energetically favored state is = 1. If at low tempera- 
tures A/b < 0.01, then A^s « 0.l7r6^, about half that at 
high temperatures. Since the corresponding value of Bq 
is then about twice that at high temperatures, h fa 0.3, 
and the quantum number for the state with the lowest 



energy becomes N — 0. 

We emphasize that Eq. (34) gives the energy of a 
ring with fluxoid number TV in the hole when there 
are no vortices in the annular region a < r < b. It 
also states that the energetically favored value of iV is 
approximately equal to Ba/BQ and hence increases as 
Bq increases. However, since the intervortex spacing 
for an infinite thin film in an applied magnetic induc- 
tion Ba is of the order of {4>q/ BaY^'^ , we expect that 
it will first become energetically favorable for a vortex 
to sit in the middle of the annular region, rather than 
to enter the hole and increase the value of N , when 
Ba = Bel ~ 'Pq/w'^1 where w ^ b — a. As pointed out in 
Ref. |3j for narrow rings [w {b ~ a)] when A/b ^ 1, 
Bel ~ (S^o/ttw^) ln(2ti;/7r^), where ^ is the coherence 
length or vortex-core radius. 



VII. LIMIT OF LARGE PENETRATION DEPTH 

In the limit A/5 ^ 1, which can be realized in small 
mesoscopic rings,SiSiS4 analytic expressions for many of 
the quantities of interest can be obtained from perturba- 
tion theory. To lowest order in the small parameter b/A, 
we may neglect the current-generated contribution Aj 
to the vector potential A in the London equation. From 
Eqs. (1) and (2) we thus obtain the sheet-current density 
J: 

From Eqs. (4) and the Biot-Savart law we then obtain 
the magnetic moment m, total current /, and magnetic 
induction at the center of the hole 6(0) generated by this 
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current: 



m 



B. 



An 



/ = 



Mo A 



In - - -^i 62 _ „2 
ZTT a 4 



277 a6 2 ^ ^ 



(36) 



(37) 



(38) 



For the case Ba — and > 0, the self-inductance 
L is obtained from Eqs. (24) and (37). Expressing L in 
terms of the same ratio as shown in Fig. 2, we obtain 



L 

Moa 



ib) (a 



2tt 



[a/b] ln(6/a)' 



(39) 



The inverse self-inductance, expressed as the dimension- 
less ratio shown in Fig. 6, is 



ai 



/ 



Mo« 
L 



2TTiA)ib) 



In 



(40) 



Expressing the magnetic field -Bi(O) at the center of the 
ring in terms of the same ratio as plotted in Fig. 4, we 
obtain 



(41) 



The magnetic moment mi expressed as the same dimen- 
sionless ratio plotted in Fig. 5, is 



6$ 



/ 



4Va 



1 - 



(42) 



For the zero-fluxoid case, we obtain the magnetic in- 
duction B2 (0) at the center of the ring by setting = 
and adding Ba to Eq. (38). Expressing the result in 
terms of the same ratio as plotted in Fig. 9, we obtain 



^2(0) ^ 1 

Ba 4VA 



1 - 



(43) 



The corresponding magnetic moment is obtained from 
Eq. (36) with <i>j = 0. Expressing the result in terms 
of the same dimensionless ratio as plotted in Fig. 11, we 
obtain 



TO2 

Too 



37r 
64 



(1) 



1 



(44) 



The total current in the ring is obtained from Eq. (37), 
again with = 0. Expressing the result in terms of the 
same dimensionless ratio as plotted in Fig. 12, we obtain 



f3i 



h 

bHa 



(l)h(i)1- 



This result is the same as that in Eq. (42), as expected 
from the identity am — Pi, proved earlier. 



For the case of flux focusing, we obtain the effective 
area Aeff = ^f/Ba by setting / = in Eq. (37). The 
result is 



A 

Aoff = ab — 
a/ 



^(&2 - g^) 
21n(Va) ■ 



(46) 



Expressing Acff in terms of the same ratio as plotted in 
Fig. 14, we obtain 



imb 



Of 

Pi 



2{a/b) \n{b/a) 
l-(a/6)2 ■ 



(47) 



The corresponding magnetic induction at the center of 
the ring B{0) is obtained by adding 6(0) [Eq. (38)] to Ba 
but using Eq. (37) and / = to eliminate $ / . Expressing 
the result in terms of the same ratio as plotted in Fig. 16, 
we obtain 



B{0) a a l/6\/ a 
Bn b ^ 6 ^ 4 VaJ V " 6 



(l-aV62) 



21n(6/a) 



(48) 



The corresponding magnetic moment to is obtained from 
Eq. (36), using Eq. (37) and / = to ehminate <&/. Ex- 
pressing the result in terms of the same ratio as plotted 
in Fig. 17, we obtain 



TO 3tt / b 
^ ~ 64 VA 



ln(6/a) 



(49) 



Finally, we calculate the quantities of interest for ring 
energies, discussed in Sec. VI. When A/6 1, the char- 
acteristic energy, cq = (l)'^/2L — 0QQ!//2^oa, becomes [see 
Eq. (40)] 



0gln(6/a) 



(50) 



which agrees with that found in Rcf. 24, allowing for the 
different system of units and definition of A used there. 
For large A, the quantity 7 = X ~ 1 = 86/3mQ://3aQ;^ — 1 
appearing in Eq. (34) becomes 



(51) 



in agreement with Ref . |2^ 

Although the results presented in this section should 
apply only when A/6 3> 1, comparison of the above ex- 
pressions with the numerical results show that they pro- 
vide good approximations already when A/6 = 1 and 
yield excellent results when A/b > 3. 



VIII. SUMMARY 

In this paper we presented a straightforward matrix- 
inversion method for the solution of the sheet current, 
vector potential, and magnetic field generated by a thin- 
film (d < A/2) superconducting ring (inner and outer 
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radii a and b) containing a trapped fluxoid in a per- 
pendicular applied magnetic induction Ba for values of 
A = X^/d ranging from zero to infinity. We used this 
method to calculate magnetic-field, current-density, and 
vector-potential profiles and numerous related physical 
quantities for three important cases: (i) <&/ > but 
Ba = (trapped fluxoid, index n = 1), (ii) Ba > 
but = (zero-fluxoid state, index n = 2), and (iii) 
$/ > and Ba > but no net current around the ring 
(flux focusing). We also calculated the Gibbs free en- 
ergy of the ring as a function of the quantum number N, 
where N = ^f/(f>Q, and the applied magnetic induction 



Ba when no vortices are present in the annular region 
between a and b. 
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